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abstract 

In  a  previous  paper  we  have  derived  a  number  of  simple  approximations  for  the 
thermodynamic  and  structural  properties  of  liquids.  Here  We  test  these  approxi¬ 
mations  by  applying  them  to  models  for  ionic  solutions  and  atomic  liquids  and 
comparing  the  results  with  those  obtained  from  Monte  Carlo  and  molecular  dynamics 
simulations.  In  particular,  the  pair-correlation  functions  for  1-1  primitive  model 
aqueous  ionic  solutions  are  calculated  using  the  EXP  approximation.  The  internal 
energy  for  2-2  salts  is  obtained  from  the  ORPA+B2  approximation.  The  pair- 
correlation  function  for  the  Lennard-Jones  fluid  is  calculated  with  the  EXP  approxi¬ 
mation.  '-In  all  cases  the  results  obtained  agree  closely  with  those  obtained  from 
computer  experiments .  Thus ,  the  ORPA+B2  approximation  for  the  free  energy  and  the 
EXP  approximation  for  the  pair  correlation  function  provide  a  theory  that  is  both 
very  accurate  and  applicable  to  a  wide  variety  of  classical  fluids.;' 
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Abstract 


In  a  previous  paper  we  have  derived  a  number  of  simple  approxi¬ 
mations  for  the  thermodynamic  and  structural  properties  of  liquids. 
Here  we  test  these  approximations  by  applying  them  to  models  for  ionic 
solutions  and  atomic  liquids  and  comparing  the  results  with  those 
obtained  from  Monte  Carlo  and  molecular  dynamics  simulations.  In 
particular,  the  pair-correlation  functions  for  1-1  primitive  model 
aqueous  ionic  solutions  are  calculated  using  the  EXP  approximation. 

The  internal  energy  for  7.-2  salts  is  obtained  from  the  ORPA  +  B2 
approximation.  The  pai r- correlation  function  for  the  Lennard- Jones 
fluid  is  calculated  wi th  the  EXP  approximation.  In  all  cases  the  re¬ 


sults  obtained  agree  closely  with  those  obtained  from  computer  ex¬ 
periments.  Thus,  the  0RPA+B2  approximation  for  the  free  energy  and 


the  EXP  approximation  for  the  pair  correlation  function  provide  «  tnv.wi  y 


that  is  both  very  accurate  and  applicable  to  a  wide  variety  of  classical 


fluids. 


I.  Introduction 


This  paper  Is  the  third  In  a  series  concerning  optimized  cluster 
expansions  for  equilibrium  classical  fluids.  The  first  article1  (I), 
presents  general  graphical  methods  and  some  simple  approximations  for 
describing  fluids  composed  of  spherical  particles.  The  second2  (II) 
generalizes  the  methods  to  apply  to  molecular  fluids.  The  object  of 
this  third  article  Is  the  study  of  the  accuracy  of  the  methods.  In 
particular,  we  apply  the  ORPA  + B2  approximation  for  the  free  energy 
[Eo.  (5.5)  of  I]  and  the  EXP  approximation  for  pair-correlation  func¬ 
tions  [Eo.  (5.7)  of  I]  to  two  different  model  fluids.  These  approxi¬ 
mations  are  the  principal  results  of  paper  I.  Their  molecular  analogs 
are  the  principal  results  of  paper  II. 

In  Section  II,  we  study  the  primitive  model  electrolytes  for  1-1 
and  2-2  aqueous  ionic  solutions.  The  pair-correlation  function  is  com¬ 
puted  for  1-1  salts  using  the  EXP  approximation.  The  internal  energy 
for  2-2  salts  is  obtained  from  the  QRPA+B2  approximation.  The  calculated 
results  agree  closely  with  those  obtained  from  Monte  Carlo  computer  ex¬ 
periments. 

The  pair-correlation  function  for  simple  liquids  is  studied  in 
Section  III.  The  EXP  approximation  Is  used  to  calculate  g(  r)  for  the 
Lenna rd- Jones  fluid.  The  results  obtained  agree  virtually  perfectly 
with  those  obtained  from  molecular  dynamics  computer  experiments. 

Some  concluding  remarks  are  made  in  Section  IV. 

Throughout  the  paper,  the  notation  and  terminology  is  the  same 


as  that  used  in  I. 
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II.  Ionic  Solutions 


In  this  section,  we  consider  the  symmetric  primitive  model  for 
aqueous  Ionic  solutions  which  represents  solvated  Ions  as  charged  hard  spheres. 
The  cations  and  anions  have  the  same  hard  sphere  diameter,  d,  and  the 
same  magnitude  of  charge,  ze,  where  e  is  the  magnitude  of  the  elec¬ 
tronic  charge.  We  let  the  subscripts  1  and  2  denote  anions  and  cations, 
respectively.  Then  the  pair  interaction3  between  an  ion  of  type  a 
and  an  Ion  of  typ^  y  Is  u^(  r)  +  u  (  r) ,  where  u^(  r)  Is  the  hard  core 
potential  (the  reference  interaction), 


ud(  r)  =  ®  t  f  <d  , 
=  0  ,  r  >  d  , 


(2.1) 


and  u  ( r)  Is  the  perturbation 

ay 

u  (  r)  =  z  z  e2/er,  r>d  .  (2.21 

CiV  u  r 

Here,  z  =  (-l)uz,  and  e  is  the  dielectric  constant  of  the  solvent, 

a 

The  number  densities,  pi  and  p2,  of  the  two  ions  are  equal  because  of 
electroneutral i ty. 

The  perturbation  potential,  u  (  r) ,  can  be  chosen  arbitrarily 

av 

for  r<d  without  changing  the  physical  properties  of  the  model  fluid. 

The  particular  functional  form  we  adopt  for  u  y(  r)  in  this  region  is 
the  one  which  satisfies  the  optimization  condition  discussed  in  I  and 
which  hence  improves  the  accuracy  of  truncated  cluster  expansions. 

Monte  Carlo  computer  experiments  have  been  performed  by  Card  and 
Valleau*  for  primitive  model  electrolytes  which  simulate  aqueous  ionic 
solutions  of  1-1  (i.e.,  z  =  1)  and  2-2  (z  =  2)  salts.  We  report  below 
results  obtained  when  the  EXP  and  OR°A  -i  B2  approximations  are  applied 
to  these  same  models. 


3- 


The  salt  concentrations  for  which  the  computer  experiments  have 
been  performed  are  sufficiently  small  that  the  Percus-Yevl ck  ( PY) 
theory  for  the  hard  sphere  reference  system  can  be  used  without  loss 
of  accuracy.  When  the  PY  theory  Is  used,  the  solution  of  the  variational 
problem  which  gives  the  optimized  u  (  r)  for  r<d  [Fg.  (4.P)  of  I] 
is  equivalent  to  the  solution  of  the  mean-spherical  model  ( MSM)  e..u«j>-,on 
for  the  fluid.5  Waisman  and  Lebowltz6  have  solved  this  problem  anr.  yt- 
ically  for  the  symmetric  primitive  model.  With  their  solution,  one  finds 
that 

u 1y(r)  -  (zaZY£2/0d)  (  2B  -  B2  r/d)  ,  r£d,  (2.3) 

and 

a..ING  "  (  12TTd3)“'  t6*  +  3x2  +  2  -  2(  1  +  2x)3''2]  .  (2.4) 

In  these  equations, 

x  =  nd  =  [8rrz2e2c/kgTe]  '^2  (2.5) 

B  =  x"2  [x2  +x  -  x(  1  +  2x)  ,/2]  ,  (2.6) 

where  c  =  pi  =pj  is  the  salt  concentration.  The  quantity  h  is  the 
familiar  Debye  inverse  screening  length. 

Using  Eqs.  (2.4)  and  (2,3),  It  Is  a  straightforward  calculation 
to  apply  the  ORPA  +  B?  approximation  [Er>.  (5.5)  of  I]  and  the  EXP 
approximation  [Eq.  (5.7)  of  I]  and  compare  the  results  with  those  ob¬ 
tained  from  Monte  Carlo  simulations*  of  ionic  solutions.  Figures  1 
and  2  show  the  pair  correlation  functions  for  1-1  salts  for  concen¬ 
trations  between  .009  and  2  M.  In  Table  I,  the  contact  values  of  the 
pair  correlation  function  are  given  for  these  solutions.  For  2-2  salts, 
the  configurational  energy  is  shown  in  Table  II.  In  all  these  cases 
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the  agreement  between  the  two  methods  Is  quite  good.  The  energy  re¬ 
sults  for  2-2  salts  are  particularly  Interesting  since  they  exhibit 
the  rapid  convergence  of  the  optimized  cluster  expansion.  Note  that 
the  perturbation  in  this  case  is  not  small.  At  hard  core  contact, 

1  uciy(  r)  )  is  greater  than  6  k^T.' 

It  should  be  noted  that  although  the  solution  of  the  MSM  equa¬ 
tions  is  used  to  obtain  the  optimized  renormalized  potential,  the  results 
of  the  ORPA  h  Bj  and  the  EXP  approximations  are  not  the  same  as  those 
of  the  MSM  equation.7  For  1-1  and  ?-?  salts,  the  MSM  predictions  for 
the  internal  energy  disagree  significantly  with  those  of  the  Monte 
Carlo  and  ORPA  +  B2  calculations  (sec  Table  II).  The  disagreement  is 
especially  striking  at  very  low  concentrations.  The  reason  for  the 
disagreement  is  that  the  MSM  Is  qual i tatlvcly  incorrect  at  low  concen¬ 
trations.  The  correct  low  density  limit  for  g^v(  r)  is  exp[  -  u  (  r) /k^T] 

for  r>d  The  MSM  equation,  however,  gives  1  -  u  (r)/k  T.  In  fact, 

—  aY  ay  B 

the  MSM  result  (for  the  1-1  and  2-2  salts  studied  in  this  section)  has 
the  unphysicol  feature  that  g  (  r)  is  negnt  ive  for  distances  at  and 
near  contact  of  the  hard  spheres  (see  Table  I).  The  EXP  approximation, 
on  the  other  hand,  has  the  correct  low  density  limit  and  predicts 
a  g  (  r)  that  is  always  non-negative.  Similarly,  the  ORPA  -t  results 
for  thermodynamic  properties  are  much  more  accurate  at  lex;  concentrations 
than  those  derived  from  the  physically  incorrect  MSM  correlation  functions. 


The  important  reason  why  the  results  of  the  optimized  cluster  theory 
are  so  accurate  is  that  the  optimization  condition  is  used  in  the  choice 
Of  perturbation  and  renormalized  potentials.  The  optimization  pro¬ 
cedure  is  much  more  significant  than  the  specific  form  of  cluster 
theory  which  we  used.  There  arc  several  ways  of  obtaining  formally 
exact  infinite  series  for  the  properties  of  fluids,  for  example  the 
Mayer  cluster  series8  and  the  mode  expans  1  on. 9  »  1  0  These  series  can  be 
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pertlelly  summed  In  a  variety  of  ways.  For  example,  chains  of  per¬ 
turbation  bonds  can  be  summed  as  In  the  Mayor  Ionic  solution  theory', 
or  generalized  chains  can  be  summed  os  In  the  y  expansion  theory. 11 
Finally,  these  partially  summed  series  can  be  reqrouped  In  a  variety 
of  ways.  For  example,  terms  of  the  same  order  In  the  concentration 
can  be  grouped  together,  or  the  y  ordering  or  T  ordering  procedures'* 
could  be  used.  Each  of  the  resulting  series  Is  formally  exact  and  so 
they  are  all  equivalent  to  each  other.  For  practical  computations, 
hcxvever,  on  Infinite  series  Is  useful  only  if  all  but  a  few 
terms  con  be  neglected  without  appreciable  loss  of  accuracy. 

Tor  series  which  contain  renormalized  potentials,  the  optimization  con¬ 
dition  Improves  the  convergence  so  well  that  for  many  problems  the 
terms  of  significance  can  be  easily  identified  and  evaluated  and  the 
others  con  be  neglected.  In  this  paper  we  have  chosen  a  scries  based 
on  the  Mayer  cluster  theory  and  a  summation  technique  which  is  suggested 
by  the  y  expansion,  but  results  of  comparable  numerical  accuracy  could 
also  be  obtained  using  a  number  of  other  types  of  series,  provided  the 
optimization  procedure  i s  emp 1  oved . 1 i 

The  advantage  of  using  an  expansion  of  the  form  used  herein  is 
that  all  the  terms  in  the  expansion  are  known  fairly  explicitly.  As 
a  result,  the  significant  terms  in  the  series  for  the  free  energy  and  g(  r)  can 

be  easily  Identified  (as  was  done  in  paper  I),  In  contrast,  the  successive 
terms  in  the  mode  expansion  and  the  y  expansion  are  not  easily  expressed 
In  ways  which  readily  allow  one  to  assess  their  importance.  Thus, 
accurate  theories  are  easier  to  obtain  with  the  cluster  expansion 
methods  used  herein  than  with  some  other  theories.  For  example  for  2-7  salts 
at  low  concentrations,  the  ORPA  + approximation  is  more  accurate  than 
the  optimized  mode  expansion10  truncated  after  the  two-mode  term, 
because  the  higher  mode  terms  contain  contributions  which  are  signifi¬ 
cant  but  which  are  difficult  to  obtain  In  simple  form. 


In  this  section  we  consider  t ho  Leonard- Jones  fluid.  The  total 
potential  energy  Is  a  sum  of  pair  potentials  of  the  6-1?  type.  Specifically, 
the  two- par tl cle  Interaction  Is 


w(  r)  =  4e  [(o/r) 1 2  -  (c/r)6]  , 


(3.1) 


where  e  mid  o  are  constants  with  the  units  of  energy  and  length,  re¬ 
spectively. 

The  analysis  of  this  fluid  begins  with  the  separation  of  the 
potential  Into  a  reference  part,  u0( r) ,  and  a  perturbation  u( r' . 

Since  It  i s  des? rable’ 5 • 1 4  for  the  reference  potential  to  be  harshly 
repulsive  (i.e,,  the  Mayer  f-function  for  u0( r)  should  closely  resemble 
the  f-function  for  a  hard  sphere  potential),  and  since  It  is  also 
desirable  for  the  perturbation  to  be  a  slowly  varying  function  of 
position,  a  very  convenient  choice  for  u0(  r)  and  u(  r)  Is 


* 

» 

uo(  0  -  w(  r)  +  e  , 

r  <  ?'/* 

0  , 

: 

«  o 

r  > 

o  , 

- 

. 

■ 

u(r)  -  -e 

r  <2'/b 

0 , 

, 

“  w(  r) 

r  >  2 6 

0  . 

(3.2) 


The  sum  u0(  r)  -1-  u(  r)  is  w(  r)  .  The  reference  potential  contains  all 
the  reoulsive  forces  in  w(  r)  and  no  other  forces.  The  perturbation  con¬ 
tains  all  the  attractive  forces.  The  motivation  for  this  particular 
separation  of  w( r)  is  discussed  in  more  detail  by  the  authors  in  pre¬ 
vious  pub i i ca t i ons.  1  5 ’ 1 4 

Since  the  reference  interaction  is  not  a  hard  core  repulsion,  the 
optimization  method  developed  in  paper  I  is  not  directly  applicable 
to  the  Leonard- Jones  system.  As  a  result,  it  is  convenient  to  introduce 
the  trial  system  in  which  the  pair  potential  is 


. TffTfP  "  \-vpw*  ^  i  . .  mm . mmfmi 
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wT( r)  =  ud(r)  •«  u(r)  .  (3.3) 

[The  subscript  "T"  Is  used  to  denote  properties  of  the  trial  system.] 

The  optimized  cluster  expansion  methods  can  be  used  to  predict  the 
properties  of  the  fluid  with  the  pair  potential  w^.(  r)  .  Once  these  prop¬ 
erties  are  known,  the  properties  of  the  Lennord-Jones  fluid  ore  straight- 
forwardly  obtained  through  the  application  of  a  type  of  cluster  expan¬ 
sion15  (which  we  have  called  the  "blip  function"  expansion),  The 
formulas  needed  to  arrive  at  the  Lennurd- Jones  properties  frem  the 
trial  properties  are  given  In  Eqs.  (11)-(14)  of  Ref.  14. 

Wo  have  used  the  EXP  approximation  for  g(  r)  [Eo.  (5.7)  of  I]  to 
calculate  gT( r) ;  and  then  with  the  blip  function  expansion  we  deter¬ 
mined  g(  r)  for  the  Lenna rd- Jones  system.  Tne  hard  sphere  radial  dis¬ 
tribution  function,  g^( r) ,  is  needed  to  perform  these  computations. 

This  function  has  been  expressed  in  convenient  analytic  form  hy  Verlet 
and  V.'eis,16  and  we  use  their  expressions.  The  optimization  condition, 

Eq.  (4.1)  of  I,  is  solved  numerically  using  the  ca 1 cula t I ona 1  procedure 
outlined  In  Section  III  of  Reference  13.  Some  of  our  results  are  graphed 
In  Figs.  3  and  4  and  tabulated  in  Tables  III  and  IV.  These  results  are 
compared  with  those  obtained  from  Verlet's  molecular  dynamics  experi¬ 
ments.17  Also  shown  for  comparison  are  the  predictions  of  the  high 
temperature  approx imati on1 3  (HTA)  and  the  optimized  random  phase 
approximation14  (ORPA). 

The  results  in  Figure  3  and  Tabic  III  describe  a  state  which  is 
near  the  critical  temperature  but  at  a  density  slightly  higher  than 
critical.  Figure  4  and  Table  IV  concern  a  high  density,  low  tempera¬ 
ture  state  near  the  triple  point,  It  can  be  seen  that  the  HTA,  ORPA, 
and  EXP  results  represent  a  sequence  of  approximations  that  converge 
very  quickly  to  the  correct  answer  and  that  the  accuracy  of  tie  EXP 
approx imat i on  is  excellent. 


T 
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IV.  Concluding  Remarks 


i 

i . 


The  results  tabulated  and  graphed  above  demonstrate  the  accuracy 
of  the  approximations  developed  In  paper  I  and  generalized  In  XI.  It 
Is  seen  that  the  0RPA  +  B2  and  EXP  approximations  are  capable  of  treat- 

i 

Ing  fluids  In' which  the  traditional  perturbation  parameters  are  large. 

While  these  approximations  are  rigorously  exact  for  low  densities  and/or 

'  I 

weak  Wnteractions^  our  results  show  that  they  rema In  accurate  when  the 

I 

density  Is  large,  -when  the  strength  of  the  perturbation  interaction  is 
large,  and  when  the  range  of  the  interaction  is  large.  Thus,  the 
ORPA  + B2  and  EXP  approximations  provide  a  comprehensive  theory  of 
classical  fluids  —  a  theory  which  is  useful  for  many  classes  of ^systems. 
Further,  the  approximations- are  computationally  simple  to  use.  The 
features  of  wide  applicability  and  simplicity  have  motivated  us  to 
cmplcvt  these  seme  methods  {along  the  lilws  of  pap^r  II)  to  study  li¬ 
quid  water.  This  investigation  is  the  subject  of  a  future  paper*. 


I 


V 
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Table  I.  Comparison  of  Radial  Distribution  Functions  at  Contact 
for  a  1-1  Symmetric  Primitive  Model  Aqueous  Electrolyte 
Solution. a 


b 

c 

^ 9+-)  MC 

( 9+-) MSM 

^9+-^EXP 

<9++>mc 

^  9++)  MSM 

(9++^EXP 

.00911 

5.22+. 40 

2.49 

4.43 

.270+. 068 

-.486 

.226 

.  1038 

3.  25+ .  07 

2.20 

3.32 

.319+.012 

-.175 

.309 

.475 

2.6?+. 04 

1.99 

2.68 

.418+.018 

.118 

.413 

1.000 

2.73+. 05 

1.90 

2.45 

.  50 5+.  27 

.358 

.522 

1.968 

2.38+. 07 

1.93 

2.45 

.706+. 028 

.631 

.669 

a.  For  these  calculations,  the  Ionic  diameter  Is  4.25  ft,  the  dielectric 
constant  of  the  solvent  Is  76.5,  and  the  temperature  Is  298°K.  Here 
g^_  Is  the  pair  correlation  function  at  contact  for  tons  of  different 

charge,  and  g  Is  the  same  quantity  for  tons  of  the  same  charge. 

TT 

MC  denotes  Monte  Carlo  results  of  Card  and  Velleau  (Ref.  4).  MSM  and 
EXP  denote  the  mean  spherical  model  and  EXP  approximations,  respectively, 

b.  The  salt  concentration  In  moles  per  liter. 

c.  Monte  Carlo  results  for  200  Ions. 


«*sa 
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Table  II.  Excess  Internal  Energy  for  a  2-2  Symmetric 
Primitive  Model  Aqueous  Electrolyte  Solution,3 

-AE/2NkBT 


b 

c 

0RPAC 

ORPA  +  Bjd 

MC6 

Af 

.0625 

1. 46 

2,02 

1.893+.017 

.13 

.250 

2.  18 

2,55 

2.473+, 019 

.08 

.5625 

2.64 

2.90 

2.822+. 008 

.08 

1.000 

2.  98 

3.  19 

3.091+.011 

.  10 

2.000 

3.38 

3.57 

3. 509+.016 

.06 

j 

£ 


i  „ 

i  ■  1 

i  j 

I  I 


a.  In  this  model,  the  tonic  diameter  is  4,25  the  dielectric 
constant  Is  78,358  and  It  is  temperature  independent,  the 
temperature  is  298. 16°K.  The  Quantity  AE  is  the  configura¬ 
tional  internal  energy  of  the  ions  and  ?N  is  the  total  number 
of  lens, 

b.  Salt  concentration  In  moles  per  liter, 

c.  Obtained  from  differentiation  of  the  ORPA  for  the  free  energy. 
For  symmetric  primitive  model  electrolytes,  the  energy  eouation 
applied  to  the  MSM  gives  identical  results, 

d.  Obtained  f rom  d i f feren t i a t i on  of  the  ORPA  +  B2  approximation 
for  the  free  energy. 

e.  Monte  Carlo  predictions  (Ref.  4). 

f.  Difference  between  the  ORPA + B2  and  Monte  Carlo  results. 
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Tab]  o  III.  Comparison  of  Radial  Distribution  Functions  for 

the  Lcnnard-Jones  Fluid  for  p*  --  .50  and  T1'  -  1.36 


r/o 

Wr,a 

90RPa(  r' 

9 EXp( 

/  \<* 
9MD<r> 

Ae 

0.  840 

0.000 

0.000 

0.000 

.000 

.000 

0.880 

0,001 

o.oo; 

0.001 

.001 

.000 

0.920 

0.053 

0.054 

0.054 

.  069 

.015 

0.960 

0.404 

0.41  1 

0.420 

.428 

.008 

1.000 

1.055 

1.087 

1.  128 

1.109 

-.  019 

1. 040 

1.579 

1.653 

1.737 

1 . 694 

-.  043 

1. 080 

1.770 

1.889 

2.007 

1.961 

-.  046 

1.120 

I  .713 

1.870 

2.005 

1.946 

-.  059 

1.  160 

1 .  :03 

1.749 

1. 868 

1.  830 

-.038 

1.200 

1.468 

1,610 

1 . 692 

1.658 

-.034 

1. 240 

1.368 

1.474 

1.521 

1.491 

-.  030 

1.280 

1.200 

1.350 

1.373 

1.335 

-.038 

1.320 

1 . 204 

1.241 

1.249 

1.227 

-.  022 

1.360 

1.139 

1.148 

1.150 

1.129 

-.021 

1.400 

1.083 

1.072 

1.071 

1.059 

-.012 

1.480 

0.996 

0.960 

0.961 

.929 

-.032 

1.560 

0.938 

0.894 

0.  898 

.  881 

-.017 

1.640 

0.904 

0.  865 

0.  869 

.  860 

-.009 

1.72C 

A  f\r\  a 

V.  <JJ  c 

0.863 

0.  £66 

.  854 

-.  01  2 

1.800 

0.899 

0.883 

0.  885 

.  894 

.009 

1. 880 

0.923 

0.920 

0.920 

.938 

.010 

1.960 

0.961 

0.972 

0.972 

.994 

.022 

2. 040 

1.007 

1.029 

1.030 

1.047 

.  017 

2.  120 

1.032 

1.060 

1.062 

1.075 

.013 

2.  200 

1.040 

1.068 

1.069 

1.073 

.004 

2.  280 

1.037 

1.056 

1.057 

1.056 

-.  001 

2.  360 

1.027 

1.037 

1. 037 

1.030 

-.007 

2.400 

1.022 

1,026 

1.027 

1.022 

-.  005 

a.  High  temperature  approximation. 

b.  Optimised  random  phase  approx imat i on 

c.  EXP  approximation. 

d.  Molecular  dynamics.  Ref.  17. 

e.  The  difference  between  the  molecular  dynamics  result  and  the  EXP 
approx  i inn  t  i  on. 
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Table  IV.  Comparison  of  Radial  Distribution  Functions  for  the 
Lennard-Jones  Fluid  for  =  .35  and  Tl!'  -  ,88 


r/a 

9HTa(  r) 

9  OR  PA' 

9EXp(  ^ 

/  \d 
W  r) 

Ae 

0.  840 

0.  000 

0.000 

0.  000 

0.  000 

.000 

0.880 

0.  000 

0.000 

0.000 

0.  003 

.003 

0.920 

0.  019 

0.018 

0.014 

0.  030 

.016 

0.960 

0.445 

0.428 

0.339 

0.412 

.  073 

1.000 

1. 813 

1.750 

1.492 

1.  51  1 

.019 

1,040 

2.918 

2.  839 

2.599 

2.  546 

-.  053 

1.080 

2.  956 

2.924 

2.  855 

2.  871 

.016 

1.  120 

2.432 

2.480 

2.  551 

2.  594 

.  04  3 

1.  160 

1.918 

2.013 

2.  108 

2.  151 

.  043 

1.200 

1.558 

1 . 648 

1.705 

1  .  744 

.  039 

1. 240 

1.290 

1.356 

1.378 

1.386 

.  008 

1.280 

1. 083 

1.119 

1.123 

1.  121 

-.  002 

1.320 

0.921 

0.931 

0.930 

0.  94  6 

.  Olo 

1.360 

0.  797 

0.788 

0.790 

0.  817 

.  027 

1.400 

0.  710 

0.686 

0.693 

0.  74  2 

.  049 

1.480 

O.  624 

0.590 

0.  603 

0.  622 

.  019 

1.560 

0.630 

0.600 

0.611 

0.  599 

-.012 

1 . 640 

0.694 

0.676 

0.681 

0.  656 

-.  025 

1.720 

0.  793 

0.786 

0.788 

0.  757 

-.031 

1.800 

0.910 

0.914 

0.913 

pom 

-.  006 

1.880 

1 . 038 

1.047 

1.047 

1 . 072 

.025 

1.960 

1.  172 

1.181 

1.  183 

1  ,  198 

.  015 

2.  040 

1.283 

1.293 

1.295 

1.  267 

-.028 

2.  120 

1.253 

1.270 

1. 274 

1 .  246 

-.028 

2.200 

1.  125 

1.  143 

1.  145 

1.151 

.006 

2.280 

0.995 

1.001 

1.001 

1.021 

.  020 

2.360 

0.901 

0.  893 

0.  894 

0.923 

.  029 

2.400 

0.  875 

0.  862 

0.  864 

0.  881 

.017 

a.  High  temperature  approximation. 

b.  Optimized  random  phase  approximation. 

c.  EXP  approximation. 

d.  Molecular  dynamics.  Ref.  17, 

e.  The  difference  between  the  molecular  dynamics  result  and  the  CXP 
approx imat ion. 
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Fl  qure  Capt  I  ons 

Figure  1.  Pal r-correlot Ion  functions  for  a  1-1  symmetric  primitive 
model  aqueous  ionic  solution.  The  Ionic  diameter  Is  4.?5  Rt  the 
dielectric  constant  of  the  solvent  is  78.5,  and  tho  temperature  is 
?98°K.  The  salt  concentration  is  .00911  h.  The  points  are  the 
results  of  Monte  Carlo  calculations  (Ref.  4)  for  POO  Ions.  The 
curves  are  the  results  of  the  EXP  approximation.  The  upper  and 
lower  sets  of  results  correspond  to  the  function  for  unlike  and  like 
charged  species,  respectively. 

Figure  ?..  Pair-correlation  functions  for  a  1-1  symmetric  primitive 
model  ooueous  ionic  solution.  Triangular  and  circular  points 
are  Monte  Carlo  results  for  .4P5  M  and  1.986  M,  respectively, 
and  the  smooth  curves  are  the  exponential  approximation  results. 

See  the  caption  of  Fig.  1  for  additional  details. 

Figure  3.  Pair-correlation  functions  for  the  Lennard-Jones  fluid  for 
a  state  near  the  critical  point.  The  open  circles  are  molecular 
dynamics  results  (Ref.  17).  The  dashed,  dotted,  and  solid  curves 
are  the  results  for  the  high  temperature  approximation,  the  optimized 
random  phase  approx imat ion,  and  the  EXP  approximation,  respectively. 

Figure  4,  Pa  I r- correla t ion  functions  for  the  Lennard-Jones  fluid  for 

a  state  near  the  triple  point.  See  caption  of  Fig.  3  for' add i t ional 


information. 
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